function rmse=FittingError(a,b,u0,T0,t,Temperature)
% 搜索数据与真实数据拟合 均方差误差
    mine_T=PDE(a,b,u0,T0,t);
    n=length(mine_T);

    %rmse = sqrt(sum((mine_T-Temperature).^2)/n); %计算均方差误差
    rmse=sum((mine_T-Temperature).^2)^0.5/n;
end

